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Abstract 

The Modified Chaplygin Gas (MCG) model belongs 
to the class of a unified models of dark energy and dark 
matter. In this paper, we have modeled MCG in the 
framework of f-essence cosmology. By constructing an 
equation connecting the MCG and the f-essence, we 
solve it to obtain explicitly the pressure and energy den- 
sity of MCG. As special cases, we obtain both positive 
and negative pressure solutions for suitable choices of 
free parameters. We also calculate the state parameter 
which describes the phantom crossing. 

Keywords Modified Chaplygin gas, K-essence, F- 
essence, G-essence. 



1 Introduction 

Several complementary cosmological observations guide 
us that our Universe is exper iencing an accel erated 



expansion in th e current era (jPerlmutter et all 11999 ; 
Riess et all 1 19981 ). From the observations of Wilkinson 
Microwave Anisotropy Probe (WMAP) satellite which 
gathered data about the cosmic microwave background 
radiation, such a cosmic acc eleration is produce d by a 
so-called dark energy (DE) (jSherwin et alll201ll ). Such 
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a new clement of the Universe, capable of accelerat- 
ing, must be, in accordance with the Friedmann equa- 
tion, have a pressure less than minus one third of the 
energy density. The notion of a cosmic medium pos- 
sessing negative pressure is not new in cosmology. The 
very first proposal for dark energy is the cosmological 
constant denoted by A. Originally proposed by Ein- 
stein to construct a theoretical model of the universe 
such that the spherical configuration of matter in the 
universe is balanced with the negative pressure of cos- 
mological constant. Thereby creating a static universe, 
in the conformity of the notions held by many scien- 
tists of that time. The origin of cosmic acceleration 
in recent time needs A but it faces serious problems of 
fine tuning (the value of A is several orders of magni- 
tude larger then estimated from the empirical results) 
and cosmic coincidence problem (the energy density of 
matter and dark energy component are approximately 
same in our presence). Then theorists looked for other 
candidates of dark energy. There has been a wide va- 
riety of theoretical models of dark energy constructed 
in the literature including quintessence, phantom en- 
ergy, Ch aplygin gas, tachyon and dilaton dark energy 
etc, see (jCopeland et all 120061) for further details. The 
quintessence and phantom energy models are based on 
spatially homogeneous and time dependent scalar fields, 
in conformity of the cosmological principle. The La- 
grangians for the quintessence (phantom energy) has 
positive (negative) kinetic energy. Where quintessence 
faces still fine tuning problem of the parameters of its 
potential, while the phantom energy gives very eso- 
teric possibilities of Big Rip and black hole evaporations 
( Jamil et aj|2008l) . 

Quintessence as a model of dark energy relies 
on the suitable choice of the potential function 
or the potential energy of scalar fields. It is also 
possible that the cosmic acceleration could appear 
due to modification of the kinetic energy of the 



scalar fields. Such modifications are termed non- 
canonical. The kinetically driven cosmic acceleration 
was originally pr oposed as a model for infl ation, 
namely k-inflation ( Armendariz-Picon et allll999i ) . and 
then as a model for da r k ene rg y, nam el y k-essence 
(Armendariz-Picon et al 



De Putter fc Linde 



Karami et a|l2 011: 
Adabi et all [20 lit 



200fi IChibal l2002t IScherrerl 12004 lYang fc Gaol 2011 : 



2001, 2000; 



Chiba et al 



20071; ICapozziello et al 



2010: 



Khodam-Mohammadi fc Taiill2010 : 
Farooq et all I2OIOI) . This model 



is free from fine-tuning and anthropic arguments. 
K-essence has been proposed as a possible means of 
explaining the coincidence problem of the Universe 
beginning to accelera te only at the present epoch 
( Malquarti et all 120031 ). Instead, k-essence is based 
on the idea of a dynamical attractor solution which 
causes it to act as a cosmological constant only at the 
onset of matter-domination. Consequently, k-essence 
overtakes the matter density and induces cosmic 
acceleration at about the present epoch. In some 
models of k-essence, the cosmic acceleration continues 
forever while in others, it conti nues for a finite duration 
( Armendariz-Picon et alluOOOl ). 

In the last years, the k-essence model has received 
much attention. It is still worth investigating in a sys- 
tematic way the possible cosmological behavior of the 
k-essence. Quite recently, a model named g-essence 
is proposed (lYerzhanov et all l2010a|) which is a more 
generalized version of k-essence. In fact, the g-essence 
contains, as particular cases, two important models: 
k-essence and f-essence. Note that f-essence is the 
fcrmionic counterpart of k-essence. 

To our knowledge, in the literature there are rela- 
tively few works on dark energy models with fermionic 
fields. However, in the recent years several approaches 
were made to explain the accelerated expansion by 
taking fermionic fiel ds as the grayitational sources of en 
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Armendariz-Picon fc Greene! 120031 )). In particular, it 
was shown that the fermionic field plays very important 
role in: (i) isotropization of initially anisotropic space- 
time; (ii) formation of singularity free cosmological 
solutions; (iii) explaining late-time acceleration. 

A very appealing proposal to describe the dark sec- 
tor are the so-called unified models. The prototype of 
such model is the Chaplygin gas. In the unified models, 



dark energy and dark matter are described by a single 
fluid, which behaves as ordinary matter in the past, 
and as a cosmological constant term in the future. In 
this sense, it interpolates the different periods of evo- 
lution of the Universe, including the present state of 
accelerated expansion. The Chaplygin gas model leads 
to very good results when confronted with the observa- 
tional data of supernova type la. Concerning the mat- 
ter power spectrum data, the statistic analysis leads 
to results competitive with the ACDM model, but the 
unified (called quartessence) scenario must be imposed 
from the beginning. It means that the only pressureless 
component is the usual baryonic one, otherwise there 
is a conflict between the constraints obtained from the 
matter power spectrum and the supernova tests. Note 
that many variations of the Chaplygin gas model have 
been proposed in the literature. One of them is the 
so-called Modified Chaplygin Gas model. It is impor- 
tant that the MCG model belongs to the class of a 
unified models of dark energy and dark matter. In this 
context, it is important to study the relation between 
MCG and the other unified m odels of dark energ y and 
dark matter. For example, in (|Mvrzakulovll201ll) rela- 
tionship between MCG and k-essence was established. 
In this paper, we have modeled MCG in the frame- 
work of f-essence cosmology. By constructing an equa- 
tion connecting the MCG and the f-essence, we solve 
it to obtain explicitly the pressure and energy density 
of MCG. As special cases, we obtain both positive and 
negative pressure solutions for suitable choices of free 
parameters. 

This paper is organized as follows. In section II, we 
introduce the F-essence formalism. In section III, we 
briefly discuss Modified Chaplygin gas and its connec- 
tion with the f-essence. In section IV, we construct a 
governing differential equation of our model and solve 
it for several special cases and a general case. Conclu- 
sion is presented in the last section. In the Appendix 
we provide the derivation of the equations of motion of 
g-essence, k-cssence and f-essence. 



2 F-essence 

Let us briefly present some basics of f-essence. Its acti on 
has the form (|Mvrzakulovll2010bl: ITsvba et afcoilbl) 



S = 



d A x, 



-g[R+2K{Y^M, 



(1) 



where ip and ip = ip^j denote the fermion field and 
its adjoint, respectively, the dagger represents com- 
plex conjugation and R is the Ricci scalar. The 
fermionic fields are treated here a s classically com- 
muting fields (see e.g. refs. (JRibas et all 12005 : 



Cai fe Wanj|2008t IWang et allboid iRakhi et all20ld 
Chimento et a]||2008HSaha et afcooi JVakili fe Sep; 



2007 



2008; Wei 2011; Dereli et al 2010; Balantekin fc Dereli 



anal 



Armendariz-Picon fc Greene! 120031 )). What 



we 

mean by a classical fermionic field is a set of four 
complex-valued space-time functions that transform 
according to the spinor representation of the Lorenz 
group. The existence of such fields is crucial in 
our work since despite the fact that fermions are de- 
scribed by quantized fermionic fields which do not have 
a classical limit, we assume such classical fields ex- 
ist and use them as matter fields coupled to grav- 
ity. A possible justification for the existence of clas- 
sical fe rmionic fields is given in the appe ndix of ref- 
erence ( Armendariz-Picon fc Greene! 12003 ). So for a 
more extensive and physically detailed discussion of 
the properties of such cla ssical fermionic fields we refer 
to th is fundamental work (|Armendariz-Picon fc Greene 
20031). Furthermore, K is the Lagrangian density of the 
fermionic field, the canonical kinetic term has the form 



y = i#™ M v> - (A^)r*# 



(2) 



Moreover, T^ = e^7 a are the generalized Dirac-Pauli 
matrices satisfying the Clifford algebra 



{ 1 », 1 V } = W V , 



(3) 



where the braces denote the anti-commutation relation, 
e^ denotes the tetrad or "vierbein" while the covariant 
derivatives are given by 

D fl tp = d fl ip-n l ^tp, D l j,i) = d l j,i) + ipn l j,. (4) 

Above, the fermionic connection £7 M is defined by 

(5) 



n, = — 9pa {r^-eMr°r s , 



with r^ denoting the Christofiel symbols. 

We work with the Friedmann-Robertson- Walker 
(FRW) spacetime given by 



ds 2 



-dt 2 + a 2 (dx 2 + dy 2 + dz 2 ), 



(0) 



For this metric, the vierbein is chosen to be 



(e^) = diag(l, 1/a, 1/a, 1/a), (e°) = diag(l,a, a, a). 

(7) 

The Dirac matrices of curved spacetime T^ are 

r° = 7 °,r J = a -y,r 5 = -i % /^r°r 1 r 2 r 3 = 7 5 , (8) 



r„ 



7 °, T j = aj j , (i = 1,2,3). 



(9) 



Hence we get 

Oo = 0, Qj = 0.5a7 J 7° 



(10) 



and 



Y = 0.5i(tpj°ip - ipj°ip). 



Finally, we note that the gamma matrices we write in 
the Dirac basis that is as 



a k \ 5 (Q I 

T = l n t I > T = I _fc n . 7 = 



(I 

" ,0 -I 



-<T k 



°J I 

(11) 



where / = diag (1,1) and the a k are Pauli matrices 
having the following form 



1 

1 



-i 

i 







(12) 



We are ready to write the equations of f-essence (In 
detail, the derivation of these equation we will present 
in the Appendix) . Here we write the final form of these 
equations. We have 

3H 2 -p = 0, (13) 

2H + 3H 2 +p = 0, (14) 

K Y ^ + Q.^HKY + k Y )ip-il K^ = 0, (15) 

KY^ + O.b^HKY + Ky^ + iK^ = 0, (16) 

p + 3H(p + p) = 0, (17) 



where 



p = YK Y -K, p = K, 



(18) 



are the energy density and pressure of the fermionic 
field. If K = Y - V, then from the system (13)-(18) we 
get the corresponding equations of the Einstcin-Dirac 
model. At last, we note that f-essence is the particu- 
lar case of g-essence for which the energy density and 
pressure are given by 



p = 2XK X + YK Y - K, p = K, 



(19) 



where X = 0.5cb 2 is the canonical kinetic term for the 
scalar field <\> (see the Appendix). 



3 Modified Chaplygin gas and f-essence 

In the cosmological context, the Chaplygin gas was 
first suggested as an alternative to quintessence and 



was demonstrated to have an increasing A (cosmo- 
logical c onstant) behavior for the evolution of the 



Hence follows that 



universe (jKamenshchik et all l2001t iBento et all 120021 ) . 
The EoS of the M CG dark energ y model was pro- 
posed by Bcnaoum (JBenaouml 120021 ) as an exotic fluid 
which could explain the cosmic accelerated expansion. 
Later on, it was shown that the EoS of MCG is valid 
from the radiation era to ACDM model (jDebnath et al 
20041 ). The MCG parameters a and B have been con- 



strain ed by the cos mic microwave background CMB 
data (iLiu fc Lill2005l) . The stable scaling solutions (at- 
tra ctor) of the Freidmann equa tion have been obtained 
in d Jamil fc Rashidlbooslliooj) . The MCG is given by 



( Benaoum 



20021) 



B 

P = A P ; 

P a 



(20) 



where A and B arc positive constants and < a. < 1. 
In the light of 3-year WMAP and the SDSS data, the 
MCG best fits the data by choosing A — —0.085 and 
a = 1.724 (|Lu et all 120081) . The dynamical attractor 
for the MCG exists at u = — 1 (where p = cop), hence 
MCG can do the phantom crossing from u > —1 to 
uj < — 1, indepen dent to the choice of initial conditions 



( Jing et alll2008l). A gen eralization of MCG was pro- 
posed in JDebnathll2007h by taking B = B(a) = B a n , 
where n and B$ are constants. The MCG is the gen- 
eralization of generali z ed Chaplygin gas p = — B / p a 
( Barreirofc^enl 120041 : ICarturan fc Finellil J2003h with 
the addition of a barotropic term. Using Eqs. (17) 
and (20), we can show th at the MCG en ergy density 
and pressure are given by (JBenaoumll2002l ) 



P 

P 



B(l + A)- 1 +Ca-^ 1+a ^ 1+A '> 
[AB(1 + Ay 1 + ACa- 3{1+a){1+A) ] 
B\B{l + A)- l +Ca-^ l+a ^ 1+A ^ 



(21) 
(22) 



where C is a constant of integration. This constant of 
integration C can be found from the condition that the 
fluid has vanishing pressure p = when a = clq : 

C = B[A(l + A)]- 1 al il+aK1+A) . 

For the modified Chaplygin gas case, the EoS parame- 
ter is 



ACa -3(i+ a) (i+A) _ B ( 1 + A yi 





B(l + A)- 1 + Car 


-3(l+a)(l+B) 


From 


(18) we get 








dK 






P 


d\nY ' 




P = 


K. 







(23) 



(24) 
(25) 



P + P 



dp dp da 



d\nY dadlnY 



so that we have 
d\uY p a 



da 

where p a 
solution 



P + P 

dp 

da ' 



(26) 



The last equation has the following 



t/ t^ f p , a da 



(27) 



where Yq is an integration constant. We can rewrite 
this expression as 



Y = Y e j '^ dC , 



(28) 



where ( = Ca~ 3 ( 1+a ^ 1+A \ Then the corresponding ex- 
pressions for the energy density and pressure take the 
form 



P = 



D + C 



[A(D + 0-B] 



D + C 



with D = B{1 + A)^ 1 . Then we have 
P + P = (l + AXiD + Q-^, 

P( - A(D + cr^ 

1 + a 
Hence we get 

d( = In C < 1 +°)< 1 + A ) (D + C) tt: 



P + P 
Finally we have 

(1 + oQA+ct „ 

Y = y ^(i+o.)<i+A) (£> + Q-t+s. 



(29) 
(30) 

(31) 
(32) 

(33) 

(34) 



4 Solvable f-essence cosmologies 

In this section, we consider some solvable f-essence 
models related with the MCG given by (20). 

4.1 The case: B = 

In this case the EoS takes the form 

P = Ap, (35) 

where p evolves like 

p = p Q a- 3(1+A \ p = const (36) 



so that 

p = A Pa a-^ l+A \ 



(37) 



On the other hand, the equations (24)-(25) for (35) 
give 



p = FY~ 



F = const. 



In this case, the pressure is 

p = K = AFY~*-. 

Comparison of (36) and (38) yields 

Y = Y Q a- 3A , Y = pf*F-Th, 
We get the scale factor as 



(38) 



(39) 



(40) 



(41) 



The behavior of Eq. (40) against redshift a -1 — 1 = z 
is plotted in figure- 1 where we see that this evolution is 
of power-law form. 

4.2 The case: A = 

Ignoring the barotropic term in MCG, we have 



V 



B 
"p* 



(42) 



It is called the genera lized Chaplygin gas (GCG) 
(JKamenshchik et all 120011 ). Recently it is proposed us- 
ing perturbative analysis and power spectrum obser- 
vational data that the MCG model is not a successful 
candidate for the cosmic medium unless A = 0, i.e. 



the usual GCG model is favored ( Fabris et all 120111 ). 



As well-known, the corresponding energy density and 
pressure are given by 



P 
P 



B + Ca-^ 1+a ^ 



■B 



B + Ca- 3( - 1+a ^ 



(43) 
(44) 



where C is a constant of integration. From (34), (43) 
and (44) we get the expressions for the energy density 
and pressure: 



P = 



B 



1-(WY)- 



p = -Bh 



l-(WY)- 



(45) 



(46) 



where W = const. The solution for Y is determined 
from (43) and (45): 



Y = W~ l C—° 



C + Ba 3(1+a) 



(47) 



The behavior of Eq. (47) is shown in figure-2, where 
we see that the kinetic energy of the f-essence increases 
and then stays constant at higher redshifts. Note that 
this conclusion depends crucially on the chosen values 
of free parameters. 

4.3 The general case 

In this section, we consider the general case when A ^ 
0, B 7^ 0. So in this case we must solve the following 

system 



p = YK Y -K, 

K = Ap-^ 
P 



(48) 
(49) 
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Fig. 1 Evolution of kinetic energy Y as a function of 
redshift z from Eq. (40). Other parameters are fixed at 
Y = 2 and A= -1. 



Fig. 2 Evolution of Y against z from Eq. (47). Other 
parameters are fixed at C = 0.5, B = 0.4, W = 0.3 and 
a = 1.2 



p = Ypy - p, 
B 

V = a p-—- 

P a 
Solving equation (48) for K, we arrive at 



p=K = EY + Y J ^dY, 



(50) 

(51) 



(52) 



where E is an integration constant. Note that if p = 
V(^,i/>) then from (52), it follows that K = EY - 
V(i/j,ip) i.e. the purely Dirac case. Eqs. (51) and (52) 
give 

which has the following solution 

(l + A)p 1+a 

where W is a constant and 
a(l + A) 



(ffy)^p«(i+") 



s = o, 



(53) 



(54) 



A + a(l + A)' 

From (55) it follows that A and a are related by 

. (n — l)a nA 

A = - -r-± — -r-^ or a 



(n — l)a + n 



(n-l)(l + A) 



(55) 



(56) 



The search of the analytical solutions of Eq. (54) is a 
tough job. So let us find some particular solutions for 
some values of n. 

4- 3.1 Example 1: n = 

It follows from (55) that this case (n = 0) realized as 
a = ot A = —1. 1) Let us first consider the case 
a = 0. Then n = and the equations (48)-(49) take 
the form 



p = YK Y - K, 
K = Ap- B 

and the equation (54) becomes 

(1 + A)p 1+a - [WY) 1 ^ -B = 0. 

Hence we write 

p= (l + Aj-^B + iWY) 1 ^}. 
and for the pressure 

P =(l + A)- 1 [-B + A{WY) 1 ^}. 



(57) 
(58) 



(59) 
(60) 
(61) 



The corresponding equation of state (EoS) parameter 
is given by 



= A- 



B(l + A) 

B+(WY)~ 



(62) 



2) Now we consider the case when A = — 1. Then 
n = and equations (48)-(49) take the form 



p = YK Y - K, 
K = -p- Bp- a . 

Hence we get 

aB\np-{l + a)- l p l+a =ln(C 1 Y)- 



(63) 
(64) 



(65) 



Consider some particular solutions of this equation. If 
a = 0, then we have 

p = ln(C 1 Yf, (66) 

and for the pressure 

p = \n(dY)- B -B. (67) 

The corresponding EoS parameter is given by 

uj = -l-[ln{CiY)}- 1 . (68) 

Second example is a — — 1 . Then for the energy density 
and pressure we obtain 

p = C 2 Yttb (C 2 = const), (69) 

and 

P=-{1 + B)C 2 Yttb. (70) 

The corresponding EoS parameter is given by 

u = -l-B. (71) 
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Fig. 3 For n = 2, the EoS parameter u»_ is plotted against 
the kinetic energy Y, while other parameters are fixed at 
B = 0.1, W = 0.2, a= 1.5 



4-3.2 Example 2: n = 1 



If n = 1 then from (56), it follows that A = 0. This case 
was considered in subsection 4.2, so we omit it here. 

4.3.3 Example 3: n = 2 

Now we consider the case when n = 2. In this case A 
and a related by 

2.4 



A = — 



or a = — - 



a + 2 1 + A' 

The equation for p (54) takes the form 

(1 + A)p 1+a - (WY) 
It has the solution 



^V(i+«) _ B = 0. 



(72) 



(73) 



(WY)-i{^- 



x l=p \fl-B(l + a) 2 (WY)- 
The pressure is given by 
a(WT)-» r 1 



}^. (74) 



Q - 



-{- 



lqp yi-B(l + af(WY) 



2(1 + °) 



. B (WYf[^- 



lT\Jl-B{l + a) 2 (WY)' 



The corresponding EoS parameter is 



r 



J - "-.™ 



CJzr = - - 



-B 



(w)-^{-L- 

L 1 + a 



lT\/l-5(l + a) 2 (^)- 



}. (76) 
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Fig. 5 For n = 0.5, the EoS parameter w_ is plotted 
against the kinetic energy Y, while other parameters are 
fixed at B = -0.1, W = 0.2, a = -1.5 



or 
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Fig. 6 For n — — 1, the EoS parameter w+ is plotted 
against the kinetic energy Y, while other parameters are 
fixed at B = -1, W = 2, a = 1.5 



Fig. 4 For n — 2, the EoS parameter w+ is plotted against 
the kinetic energy Y, while other parameters are fixed at 
B = 0.1, W = 0.2, a= 1.5 



Iii figures (3) and (4), we have plotted the EoS param- 
eter and it is shown that subnegativc values of u arc 
permissible in our model. This corresponds to f-essence 
MCG behaving as pha ntom energy which c auses super- 
accelerated expansion (jDebnath et all 120041 ) . 



4.3.4 Example 4-' n = 0.5 

If n = 0.5 then A and a satisfy the relation 
a A 



A 



or a 



1 + A 



a-1 

Eq. (54) becomes 

(1 + A)p 1+a - (WY)^ p°- 5{1+ ^ -B = 
This equation has the following solution 

1-ar 



(77) 



(78) 



(WY)i{- 



l±il 



45 (WY)- 1 -^ 



1-a 



r 



(79) 



The pressure is given by 
a(WY)i (1-a 



P 



n 



1±a/1 



45 (WY)-^ 



-B 



(wYy 1 ^ 



1 — a 

1 — Q 



2 

r 



i±«/i 



\B 



(WY)- 



1 — a 
The corresponding EoS parameter reads 

a (WY)^r 



2a 



(80) 



a — 1 



1-a 



1T\ 1 



45 -(Rry)-^ 



1 — a 



(81) 



-0.8119 - 
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Fig. 7 For n = — 1, the EoS parameter w+ is plotted 
against the kinetic energy Y, while other parameters are 
fixed at B = -3, W = 0.8, a = -1.5 



In figure-5, we have plotted the EoS parameter against 
the kinetic term. Here we chose a < which corre- 
sponds to the polytropic term added in the barotropic 
equation of state. Such a f-essence polytropic EoS can 
also cause the super-acceleration. 



4-3.5 Example 5: n = 

Our next example is n 
relation 



.4 



2o 



2a- 1 
Eq. (54) becomes 

(l + A)p 1+a 

It has the solution 



2(1 + A) 



(M/y)-^ P - (1+Q ) 



■1. Then A and a satisfy the 

(82) 



B = 0. 



P = 



B± J B 2 + 4{1 + A)(WY)- 



2(1 + A) 
The pressure is given by 



(83) 



(84) 



P 



A 



B± J B 2 + 4(1 + A)(WY)- 



2(1 + A) 



B 



rB±JB 2 +4(l + A)(WY)~ 



2(1 + A) 

The corresponding EoS parameter is given by 
2a B 



"(85) 



a; q 



2a -1 

B + slB 2 



-(wyy 



1 



2a -1 



(WY)~ 



(86) 



In figures 6 and 7, we plotted the above state parameter 
against kinetic energy. We choose positive and negative 
values of a for the sake of completeness. It is apparent 
that the state parameter achieves subnegativc values 
showing the viability of our dark energy model. 



5 Conclusion 

In summary, we have modeled modified Chaplygin gas 
in f-essence cosmology. The use of MCG is useful as 
a tool of explaining dark energy and dark matter in 
a unified manner, while f-essence cosmology essentially 
suitable to describe cosmic acceleration only at present 
time. Thus the correspondence of MCG with f-essence 
is useful in learning how these two models are connected 
to each other. We studied this link by constructing a 



differential equation (54) connecting the MCG and the 
f-essence. We solved it to obtain explicitly the pres- 
sure and energy density of MCG. We observed that f- 
essence MCG has one additional free parameter namely 
W along with A, B, a. As special cases, we obtain 
both positive (37), (80) and negative (46), (68), (72) 
pressure solutions for suitable choices of free parame- 
ters. The negative pressure solution is essentially useful 
for cosmic expansion with acceleration. Prior to this, 
we studied the model with barotropic and generalized 
Chaplygin gas equation of states. 

The present work is concerned only with the cor- 
respondence of MCG with the f-essence. However it 
would be more interesting to study the dynamical fea- 
tures of this model. For instance, it is interesting to 
check whether such a model is useful to model infla- 
tion at earlier cosmic epoch and cosmic acceleration at 
the present time. Next it will be important to compare 
this model with the observational data and constrain 
its free parameters, particularly W. Also observational 
constraints on the state parameter ui will indicate if it 
corresponds to cosmological constant, quintessence or 
phantom energy. We finish our work in the hope that 
merits and demerits of f-essence will be clear in the fu- 
ture investigations. 

We note that this paper is the logical continuation 
of (|Mvrzakulovl 120111) . where the relation between k- 
essence and MCG was studied. The action of k-essence 
reads as 



S 



d i xy/^j[R + 2K(X, (/>)], 



(87) 



where the kinetic term X for the scalar field (f> (for the 
FRW metric) reads as 



X = 0.. 



(88) 



The corresponding equations for the FRW metric look 
like 



3H 1 - p = 0, 

r2 



(89) 
2H + 3H'+p 0, (90) 

K x 4> + (k x + 3HK x )4>-K^ = 0, (91) 

P + 3H(p + p) = 0, (92) 

where the energy density and pressure of the scalar field 
is given by 



p = 2XK X -K, p = K. 



(93) 



In this case, the common system of equations for the 
k-essence and MCG has the form 



p = 2Xp x - p, 
B 

V = a p-—- 

p a 



(94) 

(95) 



The compat ibility condition fo r the equations (|94]) - (|95 
is given by (|Mvrzakulovll20 1 II ) 



(1 + A)p 1+a - (WI)^p n(1+a) 



B = 0, 



(96) 



where W = W{<f>) and 

a(l + A) 

T) = 

A + a(l+A)- 
In ( MvrzakulowOl II ) . different type solvable k-essence 



(97) 



cosmologies compatible with the MCG model are found 
for the different values of n. 
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6 APPENDIX. The derivation of the equations 
of motion of g-essence, k-essence and 
f-essence 

In this Appendix we would like to present the deriva- 
tion of the equations of motion for the f-essence action 
(2.1) that is the system (13)-(17). But, as f-essence is 
the exact particular case of g-essence, we consider more 
general case and give the derivation of the equations of 
motion for g-essence. Let us consider the following ac- 
tion of g-essence 



S = 



d A x, 



-g[R + 2K(X,Y,4>^,ij)}, 



(98) 



where R is the scalar curvature, X is the kinetic term 
for the scalar field <p, Y is the kinetic term for the 
fermionic field tp and K is some function (Lagrangian) 
of its arguments. In the case of the FRW metric 



ds 2 = -dt 2 + a 2 (dx 2 + dy 2 + dz 2 ), 
R, X and Y have the form 



, a a 
R=6 -+^ 



X = O.50 2 



Y = 0.5i(V>7V -VW), 



(99) 



(100) 



(101) 



(102) 
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respectively. Substituting (100)-(102) into (98) and in- 
tegrating over the spatial dimensions, we are led to an 
effective Lagrangian in the mini-superspace {a, (p, ip, ip} 



L = -2(3ad 2 - a 3 K] 



(103) 



Variation of Lagrangian (103) with respect to a yields 
the equation of motion of the scale factor 



2aa + a 2 + a 2 K = 0. 



Now by varying the above Lagrangian (103) with re- 
spect to the scalar field <p we obtain its equation of 
motion as 



Ky 



K x c 



3-K x ^-K 4> = 0. 
a 



(105) 



At least, the variation of Lagrangian (103) with respect 
to ip, ip that is the corresponding Euler-Lagrangian 
equations for the fermionic fields give 

K Y j°ip + l.h-K Y fiP + 0.5X Y7 > - iKj = 0, (106) 



(117) 

Also we note that g-essence admits two particular cases 
(reductions): k-essence and f-essence. In fact, let ip = 0. 
Then Y = and the system (110)-(115) becomes 



(104) K x 4> + (Ax 



3H 2 - p = 


= o, 


(118) 


2H + 3H 2 +p = 


= o, 


(119) 


IHK X )4>-K4 -- 


= o, 


(120) 


p + 3H(p + P ) -- 


= o, 


(121) 



where 



p = 2XK X - K, p = K. 



(122) 



It is the system of equations of k-essence. Now let 
us consider the case when (f> = that is the purely 
fermionic case. Then X = 0, K = K (Y, ■ijj , ip) and the 
system (110)-(115) takes the form of the equations of 
f-essence that is (13)-(17). 



K Y ipi a + 1.5-K Y ifa + 0.5AV^7° + iKj, = 0. (107) 
Also, we have the "zero-energy" condition given by 

Lah + L^ + L^ + L^-L = {) (108) 

which yields the constraint equation 

YK Y - K = 0. (109) 



-3a- 2 a 2 + 2XK X 



Collecting all derived equations (104) - (107) and (109) 
and rewriting them using the Hubble parameter H = 
(lna) t , we come to the following closed system of equa- 
tions of g-essence (for the FRW metric case): 

3H 2 -p = 

2H + 3H 2 +p = 

Kxt+ikx + ZHKxti-Kj, = 

K Y iP + 0.5{3HK Y + K y )iJj -i-/°K^ = 

K Y {p + 0.5(3HK Y + k Y )tP + iK,j,j = 
p + 3H(p + p) = 

Here 

p = 2XK X + YK Y - K, p = K 



0. 


(110) 


0, 


(111) 


0, 


(112) 


0, 


(113) 


0, 


(114) 


0. 


(115) 



(116) 



are the energy density and pressure of g-essence. It 
is clear that these expressions for the energy density 
and pressure represent the components of the energy- 
momentum tensor of g-essence: 



Too = 2AA' ; 



YK Y - K, Tn = T 22 = T 33 = -K. 
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